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ABSTRACT 

Aims. We aim to unveil the most massive central cluster hlack holes in the universe. 

Methods. We present a new search strategy which is based on a black hole mass gain sensitive ’calorimeter’ and which links the 
innermost stellar density profile of a galaxy to the adiabatic growth of its central SMBH. In a first step we convert observationally 
inferred feedback powers into SMBH growth rates by using reasonable energy conversion efficiency parameters, e. In the main part 
of this paper we use these black hole growth rates, sorted in logarithmically increasing steps encompassing our whole parameter 
space, to conduct A-Body computations of brightest cluster galaxies with the newly developed Muesli software. For the initial setup 
of galaxies we use core-Sersic models in order to account for SMBH scouring. 

Results. We find that adiabatically driven core re-growth is signihcant at the highest accretion rates. As a result, the most massive 
black holes should be located in BCGs with less pronounced cores when compared to the predictions of empirical scaling relations 
which are usually calibrated in less extreme environments. For efficiency parameters e < 0.1, BCGs in the most massive, relaxed and 
X-ray luminous galaxy clusters might even develop steeply rising density cusps. Finally, we discuss several promising candidates for 
follow up investigations, among them the nuclear black hole in the Phoenix cluster. Based on our results, it might have a mass of the 
order of 10*'Mo. 

Key words, galaxies: elliptical and lenticular, cD - quasars: supermassive black holes - methods: analytical & numerical 


1. Introduction 

The detection of supermassive black holes (SMBHs) at the high- 
mass end of the mass scale is of tremendous importance for con¬ 
straining limitations of empirical scaling relations and thus for 
the understanding of the related evolution and growth of both 
SMBHs and cosmic structures like galaxies, galaxy groups or 
even clusters of galaxies. 

SMBH scaling relations link different properties of galax¬ 
ies such as total luminosity and bulge mass (|Kormendy & Rich-| 
stone||1995[ [McLure & Dunlop||2002[ Marconi & Hunt 2003 [ 


Haring & Rix||2004f McConnell & Ma 201 3|l, the velocity dis 


persion, cr, (|Ferrarese & Merritt|2000 Gebhardt et al.|20d0{|Gul-| 

tekin et al.|2009| [McConnell & Ma||2013| l, and the central sur¬ 

face brightness profile (Lauer et al.|2007a [Kormendy & Bender 


2009 1 [Rush et al.|2013 |Dullo & Graham|2014] l to the mass of 

the nuclear SMBH. Currently, there is a debate whether the em¬ 
pirical M, - cr relation is a consequence of black hole feedback 
regulated galaxy growth (e.g. |Silk & Rees|19Mt|Fabian|2012 1 or 
results from stochastic averaging over numerous galaxy coa es- 
cences in a cosmological merger tree ( Jahnke & Macci5|201 1[ ). 

The action of binary or multiple black holes after merger 
events of their host galaxies is suspected to be responsible for 
the formation of shallow central brightness profiles of luminous 
galaxies (Ebisuzaki et al.|[l991[ Milosavljevic & Merritt||2001 


Merritt 


2006|l with V-band magnitudes My < -21 (Lauer et] 


al.||2007b I. This scenario is supported by [Thomas et al. ( 2014|l 


who found a tangentially biased velocity distribution within the 
break radius, Rbreak, which is located at the maximal curvature of 
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the surface brightness profile. In the picture of SMBH merging, 
^break forms out, along with a tangentially biased velocity dis¬ 
tribution, by the preferential scattering/removal of stars on ec¬ 
centric orbits with low angular momentum. Therefore, the tight 
relationship as indicated by scaling relations between M, and 
properties of the central galactic brightness profile seems rea¬ 
sonable. 

In a few particular and extreme cases like MS0735 and 
Abell 85 ( McNamara et al.|2009 Lopez-Cruz et al.|2014[ ), scal¬ 
ing relations which are usually calibrated in local and “ordi¬ 
nary” galaxies yield widely different black hole mass predictions 
by factors of 10 - 100. By using the M, - Rbreak (|Rusli et'af 


2013|l, M, - Rcusp (Lauer et al. 2007a > and M, - Ldef (Kormendy 


& Bender 2009 1 relations, where Rcusp is the cusp radius of a 
Nuker profile and Ldef is the total luminosity difference between 
measured and extrapolated light profile, McNamara et al. ( |2009| l 
and |L6pez-Cruz et al.| ( |2014| l found indications for extraordinary 
massive SMBHs. The central black hole masses in MS0735 and 
Abell 85 might be as high as 10"Mo as long as their shallow 
density profiles are carved out by black hole merging processed 
The estimates of these three scaling relations are significantly 
above the predictions of the M, - cr and M, - L relations. There 
is tentative evidence that it is the M, - cr (and M, - L) relation 
which loses its predictive power (throug h a vertical tilt) above 
cr = 270 km s ' (Kormendy & Ho 20131. 

Unrelated studies give also rise to skepticism concerning the 
validity of the M, - cr and M, - L relation when dealing with 


' At least for the BCG in Abell 85, Bonfini et al. ( 20 15^ show that the 
central brightness profile can be approximated by a Sersic model with a 
small index n and without a depleted core. 
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the most extreme central cluster galaxies: By (i) using the fun¬ 
damental plane of accreting black holes, Hlavacek-Larrondo et 
al. ( 2012| l found evidence for the presence of ultramassive black 
holes (M, > IO^^Mq) being ^ 10 times more massive than pre¬ 
dictions based on the M, - cr and M, - L relations. Furthermore 
(ii), the case for exceptionally massive black holes is strength¬ 
ened by active galactic nuclei (AGN) with quasar like powers in 
form of mechanical feedback but without any detectable X-ray 
emission. The central active black holes in these cool core galaxy 
clusters must operate at highly sub-Eddington accretion rates 
to avoid detection and are most likely ultramassive (Hlavacek- 
[Larrondo & Fabian|201 1| |. 

Central cluster galaxies which are located in the most mas¬ 
sive (cool core) galaxy clusters are therefore promising candi¬ 
dates for hosting the largest black holes in the universe. In¬ 
deed, stellar dynamical mass measurements of even less extreme 
BCGs in the nearby Com^ and Leo galaxy cluster as well as 
lenticular galaxies unveiled black holes around M, - IO'^Mq 
( [M cConnell et al.|2011[[van den Bosch et al.|2012| l. 

This paper investigates the impact of hot cluster gas accre¬ 
tion on the dynamics and density profiles of BCGs. We relate 
the old idea of adiabatically growing black holes ( |Young|1980| l 
in combination with core depletion by merging SMBHs to create 
a black hole mass gain sensitive “calorimeter”. This tool offers 
a way to test SMBH accretion models, improve the understand¬ 
ing of scaling relations which depend on the surface brightness 
profile and help to unveil the most massive cluster black holes. 
If combined with the total feedback energy of a galaxy clus¬ 
ter, the inner slope and break radius, f?break, of its central clus¬ 
ter galaxy and SMBH mass, the “calorimeter” might be used to 
uniquely determine the SMBH mass gain during the early quasar 
era, through subsequent merging activity and the accretion of hot 
cluster gas. Furthermore, it can be used to infer the time averaged 
binding energy conversion efficiency parameter, e. 

In Sectionj^we first compile observationally motivated black 
hole growth rates. These rates also cover the extremes, from 
a mere spin powered feedback scenario (dM,/df a; 0) up to 
SMBH growth with small e in X-ray luminou s galaxy clusters 


(Lx = 10"^^ - 10‘^^ erg s ') like the Phoenix- (McDonald et al. 


[2012} or RX J1347.5-1145 cluster ( jGhti & Schindler|2004| ). In 
the main part of this paper (Section we perform computa- 
tions with the newly developed Muesli software ( jBrockamp et 
|al.|2014| l in order to evaluate the impact of the (compiled) SMBH 
growth rates on the innermost surface brightness profiles. In Sec- 
tionj^we discuss the obtained results, present a few cluster can¬ 
didates which are recommended for follow-up observations and 
describe future applications and improvements of the “calorime¬ 
ter”. Our main findings are summarized in Section]^ 


2. Accretion models 

In this section we present central cluster black hole growth mod¬ 
els which serve as the initial setup for our A-Body computations 
in Section]^ They are obtained from the observed AGN feed¬ 
back power in combination with reasonable accretion efficiency 
parameters (Section HD- Four mutually independent methods 
are reviewed in ascending order with regard to their mass accre¬ 
tion rate. Due to the vast amount of literature our compiled list 
is far from being complete, but it covers the extreme scenarios. 


A consistency check is made in Section 2.3 by comparing the 


observationally motivated black hole growth rates with classi¬ 
cal Bondi models in order to obtain reasonable initial black hole 
masses for our A-Body computations. Readers which are already 
familiar with deposited AGN feedback powers can directly skip 
to the new results presented in Section]^ 

2.1. Conversion of feedback power into the SMBH growth 
rate 

The deposited AGN feedback (jet) power, Pagn, is converted 
into cluster black hole growth rates, M,, by using the relation: 


dM. . 1 - e 

M, - —7— - (1 - e)Macc - —^Fagn 


df 


ec^ 


( 1 ) 


Here, Macc = Fagn^ c""' is the accretion rate and e is the bind¬ 
ing energy conversion efficiency parameter. In order to account 
for uncertainties related to e, we use the six efficiency param¬ 
eters, e - 1.0/0.42/0.2/0.1/0.057/0.01, which span over two 
orders of magnitude. They are motivated as follows: 

- e = 1.0/0.1: These values are obtained from the elec tromag- 
netic jet efficiency formula, e =» 0.002 ■ (1 - |a|)^* (Hawley 


& Krolik 20061, in combination with black hole spin val¬ 


ues a = 0.998 and a - 0.98. While a = 0.998 is the the¬ 
oretical limit which can be obtained in accretion processes 
( Thorne|[l974 l, a = 0.98 corresponds to the averaged spin 
of 48 luminous A GN at high redshift (z > 1.5) with masses 
M, > 3 ■ IO^Mq ( Trakhtenbrot||2014 1 . Their dimensionless 
spin parameter, a - 0.98, is derived from deduced radia¬ 
tive efficiency parameters of individual accretion disks. Ad¬ 
ditionally, e = 0.1 is often used as the canonical value for the 
energy conversion efficiency (e.g. [Frank et al.|1992||. 
e = 0.42/0.2/0.057: The value e = 0.057 (e = 0.42) is 
the maximal conversion efficiency of a (co-rotating) circu¬ 
lar accretion flow (within the equatorial plane) around a non 
spinning (maximally spinning) black hole (e.g. [Novikov'S] 
Frolov|1989|l. The parameter e = 0.2 lies in between. 
e = 0.01: The kinetic jet/feedback efficiency parameter de¬ 
rived from AGN observations ([Kording et al.|2008| [Merloni 
[& Heinz|2008t . 


^ We note that the Coma cluster does neither host a cool core 
( [Simionescu et ar|[2013[ l nor is it dominated by one central clus¬ 
ter galaxy. Instead, there are three dominating galaxies which are 
NGC 4889, NGC 4 874 and NGC 4839. T hey will likely coalesce some¬ 
time in the future ( [Gerhard et al.[[2D07[ l. Their merger will give 


we used direct SMBH mass estimates from|Magorrian et al.[([1998J for 
NGC 4874, [McConnell et al.[ ( |20TT] ( for NGC 4889 and obtained a mass 
of M. 3! 6.5 • IO^Mq for NGC 4839 with the help of the M. - Rbieak re¬ 
lation from [Rusli et al. j2013^ and a graphically evaluated break radius 
from Figure 3 i n [Jordan et al.[j2004). Ma ss loss through gravitational 
wave emission ( [Tichy & Marronetti|2008^ but also stellar accretion in 
non-spherical potentials i Liu & Chen 2013 Vasiliev 2014[( is neglected. 


2.2. Growth rate based on observationai modeis 

By using the empirical relationship between mechanical work 
and jet luminosity at 1.4 GHz from [Cavagnolo et al.[ ( [2010| l, [Ma| 
et al. ( 2013| l obtained a time and galaxy cluster averaged AGN 
power of Pagn = 3 ■ 10^"^ erg s '. The corresponding SMBH 
growth rate is: 


M.iow = 0.005 


Moyr 


( 2 ) 


Equation defines our low mass accretion model as it is 
strongly exceeded by AGN in some of the most massive cool 
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core galaxy clusters which are the most promising candidates 
for hosting the heaviest black holes. 

X-ray observations of these clusters demonstrate that the de¬ 
posited AGN feedback power required to inflate X-ray cavities 
within the intra cluster gas (ICM) exceeds Pagn = 10"^^ erg 
Hlavacek-Larrondo et al.| ( 2012| l report cavity powers being as 
high as Pagn = (4 - 8) ■ lO’^erg s ' for individual cases of 
massive galaxy clusters. Due to a duty cycle (i.e. fraction of 
time where the central cluster black hole is active) of at least 
60% - 70% ( |Dunn & Fabian|[2006t [Birzan et al.||2012] ), we as¬ 
sume a conservative value of Pagn = 3 ■ 10"*^ erg s“^ which is 
also found in (Hlavacek-Larrondo et al. 2015| l. This corresponds 
to 


M.med = 0.05 


1 


Moyr 


(3) 



Medium: Pagn“^‘^® 


which defines our medium mass accretion model. The assump¬ 
tion of linear mass increase is justified by cavity energet 
ics which seemingly do not evolve with redshift ( [Hlavacek 
[Larrondo et al.|2012| l. 

But even the value defined in Equation ([^ has to be regarded 
as a lower limit. Neither is the feedback energy fully deposited 
inside the cool core nor are cavities the only source of ICM heat¬ 
ing. Energy can be transferred out to much larger radii via large 
scale shocks or sound waves. One spectacular example is the 
aforementioned system MS0735 with its bipolar 200 - 240 kpc 
X-ray cavities surrounded by Mach 1.26 shocks waves ( |Van-| 
|tyghem et al.||2014) l. They extend to at least twice the size of 
the cooling radius. And indeed, studies which calculate the total 
injected AGN feedback energies out to large radii o btain time 
averaged jet powers as high as Pagn = 10^^ erg s“' (Mathews 
|& Guo|201 lUChaudhuri et al.|2012| l. These values are deduced 
from comparing observed potential energies of virialized mas¬ 
sive clusters (Myir ~ IO'^Mq) to those obtained from simula¬ 


tions without feedback. Eollowing Mathews & Guo (20111, the 


energetic demands and high accretion rates of these central clus¬ 
ter black holes can not be lowered by splitting them into mul¬ 
tiple sources. In most cases. X-ray cavities and shocks are as¬ 
sociated with the central BCG and not its satellites. Significant 
pre-heating of baryons prior to galaxy cluster assembly is also 
disregarded b y|Mathews & Gu^(|2011| l. 

Eollowing Churazov et al. ( |2002| l, similar (or even higher) 
AGN feedback powers are also required for balancing the self- 
cooling of the gas in the most X-ray luminous cool-core clus¬ 
ters approaching Lx ~ 10"^^ erg s“* like MACSJ1447.4-1-0827 
dHlavacek-Larrondo et al.]|20T2l l, RX J1347.5-1145 dGitti &| 
|Schindler|2004| l or the Phoenix cluster dMcDonald et al.|2012| l. 
Integrated over galaxy cluster lifetimes of seven to eight billion 
years and by assuming various efficiency parameters (from as 
high as unity for black hole spin powered models down to a few 
percent for Schwarzschild black holes, see Section [ZT] i this re¬ 
sults in up to AM, IO'^Mq of accreted mass. Eeedback ener- 

erg s * define our high mass accre- 



Te [Gyr] 


gies as high as Eagn = 10' 
tion scenario: 


,46 


M.Mgh = 0.18 


l-e 


Moyr 


(4) 


Central cluster black holes growing at such prodigious rates 
might pose a challenge for theories which assume that significant 
SMBH growth ceased early on in accordance with cosmic down¬ 
sizing of the AGN activity ( Hasinger et al.|2005] Kelly & Shen 
2013^ . Alternatively, the feedback might be powered mainly by 
trapped spin energy (e ^ 1) or central cluster black holes might 


Fig. 1: The SMBH growth rate is obtained from the central cluster AGN 
feedback power which is constrained by means of different strategies 
and cluster samples. The colors and line types indicate the binding en¬ 
ergy conversion efficiency parameters which must be adopted for the 
conversion (Equation[TJ. The cosmic time is defined to be (z = 0) = 0 
and we assume that feedback/accretion started operating at redshift 
z = 1 (8 billion years ago). SMBH masses can be obtained from the 
relation M. (t) = M.init + AM. (t). 


have reached dimensions where even AM, > IO^Mq of accreted 
gas represents only a tiny fraction of their total mass ( |Hlavace^ 
ILarrondo et al.|[2012| l. We note that the latter idea is in accor¬ 
dance with the existence of individual cases of extremely mas¬ 
sive AGN in the early universe (Ghisellini et al.|2009j Trakhten^ 


brot|2014[ [Banerji et al.|[20T5| |Wu et al.||2015[]jun et al.||2015| l 

and is conform with the upper limit of the quasar/black hole 
mass function (BHME) ( [Kelly & Shen|2013l their Chapter 3.3). 
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The obtained BHMF is calibrated by the virial mass estimates of 
58.000 Type 1 quasars. 

The observation-based models are graphically illustrated in 
Figure [T] Note that models with e = 1 (AM, = 0) are not plotted 
subject to the logarithmic spacing of the y-axis. With regard to 
the enormous spread in e (Section [2T]i, the accumulated SMBH 
mass over 8 Gyr (i.e. the typical lifetime of galaxy clusters), 
varies from AM, = 0 up to several IO'^Mq. Consequently, the 
growth of the SMBH can not be uniquely constrained from the 
deposited feedback energy. Therefore, it is important to under¬ 
stand the influence of the growing black hole on its immediate 
surroundings. This leads us to the idea of the black hole mass 
gain sensitive “calorimeter" which is part of Section]^ 


2.3. Consistency check: Classical Bondi models 


Through a comparison with theoretically justified models we can 
explore whether the initial black hole masses which are chosen 
for our N-Body computations and which can not be directly de¬ 
duced from the observed feedback powers are realistic. We focus 
on galaxy clusters with relaxed gas atmospheres and ignore cold 
gas accretion scenarios subject to the absence of a positive cor¬ 
relation between jet power and molecular gas mass ( |McNamara| 
|et al.|[20TT] l. We assume that accretion flows are not disturbed 
through shock heating by collisions and consider only well es¬ 
tablished models, i.e. classical Bondi flows ( |Bondi||T952| l with 
some improvements from more recent investigations. The usage 
of sophisticated (analytical) models like radiating Bondi flows 
( [Mathews & Guo|2012| l is left for future studies. 

SMBH growth rates related to classical Bondi accretion 
models are shown in Figure We do not consider the influ¬ 
ence of the growing black hole on the accretion rate but discuss 


some consequences this might have in Section 4.3 Our setup in¬ 
cludes three initial black hole masses with M, = 1,2,4 -IO'^Mq. 
The highest value corresponds to the upper quasar mass obtained 
from the black hole mass function ( [Kelly & Shen| 2013 i. Indi¬ 
vidual cases of black holes with such masses at high redshift are 
also found in [Trakhtenbrot[ pO 14| l and [Jun et al.[ ( [20l3| ) and are 
required to successfully model (but not strictly unambiguously) 
the spectral energy distributi on and total luminosity of the blazar 
S5 0014-H813 at z = 3.37 ( [Ghisellini et'Z]|2009|p 1 The three 


initial SMBH masses are discriminated by using different line- 
types in Figure]^ 

Apart from some modifications which are specified in the 
text, we make use of the same strategy as [Allen et'af p006 1 
and use their equations and parameter setting for calculating the 
black hole growth rate through Bondi accretion: 

M.b = 'I (1 - e) 47tA (GM.f (5) 

Here, rj - 0.2 is a conservative correction parameter for rotat¬ 
ing and viscous gas atmospheres ( [Narayan & Fabian|201 1| |. The 
term (1 - e) = 0.9 guarantees energy conservation (assuming 
e - 0.1) and the coefficient A - 0.25 is related to the adiabatic 
index 7 = 5/3. The sound velocity at the characteristic Bondi or 

accretion radius, = 2GM,/c^, is Cj = yjyk^rjjpm^. The 
parameter kg is the Boltzmann constant, T is the gas tempera¬ 
ture, p - 0.62 is the average atomic weight and irip is the mass 
of the proton. The gas density, p = 1.13«<,mp, which is also taken 
from [Allen et al.[ ( |2006[ ) can be calculated from the measurable 
electron density, n^. 


^ Note that blazars are highly beamed sources, so for every blazar there 
should be several non-beamed AGN with similar properties ([Ghisellimj 


et al. 2009 1 . 
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Fig. 2: Black hole growth models related to the classical Bondi accre¬ 
tion scenario. Although, the rates differ widely and predict extreme up¬ 
per values, they lie within the parameter range of the observation-based 
models which are plotted in Figure[^ In case of very massive seed black 
holes around M, = IO'^Mq at z = 1, Bondi accretion alone is sufficient 
to explain the deposited feedback energies even in the most extreme 
cool core clusters. 


The remaining free parameters, rig and T, which are required 
for our growth models are taken from a cluster sample average. 
To be more precise, we calculate the arithmetic mean values 
of He and T at Ra from Table 4 in [Russell et al. ( 2013 1. Their 
sample includes the same galaxy clusters as Allen et ak] ( [2006| l 
plus four additional objects. Two sets of parameters (distin¬ 
guished by different colors in Figure]^ are chosen. The first one. 
He — 0.52cm“^, T — 4.9 ■ 10®K, represents the mean, whereas 
the second one, = 0.26cm~^, T = 5.1 ■ 10®K is obtained from 
the upper (T + AT) and lower («<, - Arig) Icr mean values. In this 
way, we tend to underestimate the true SMBH growth rate. 

Uncertainties related to Ra (which depends on the assumed 
M,) and the parameters tig (Ra) and T (Rg) are compensated by 
a constant ratio cj^p oc PF 5/2 _ ^onst (in Equation!^ which is 
independent of radius, r, for y = 5/3, P = nksT and - 

const. It is only important that rig and T are measured deep within 
the potential of the central host galaxy because of its influence 
on the Bondi accretion flow. 

Evidently, all analytical-based models which are presented 
in Eigure lie within the parameter range of the observation- 
based models (Section [2. 2[ i. This indicates that Bondi accretion 
alone is sufficient for powering the most extreme central cluster 
AGN hosting the largest SMBHs. Eurthermore, it shows that ini¬ 
tial black hole masses as high as M, = (1 - 4) ■ IO'^Mq at the 
onset of hot cluster gas accretion are reasonable. 


3. Cusp formation due to SMBH growth 

Initial conditions and results of A-Body computations are pre¬ 
sented in this section. The idea is to numerically investigate, 
for a sample of representative black hole growth models (Sec¬ 
tion]^, the adiabatic response of the red/old stellar density pro¬ 
file of the host galaxy to the mass which becomes trapped by the 
black hole. Under the restriction that the duration of the SMBH 
growth phase is larger than the orbital timescale, the term ’adia¬ 
batic’ refers to the increase of the central density while simulta- 
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neously the angular momentum distribution remains conserved. 
This idea goes back to Young ( |1980| l and was later used by |van| 
der Marel ( 1999|l to explain and predict the observed power law 


cusps of early type galaxies and central SMBH masses, [van der| 
Marel ( \999\ assumed that initial galaxies already started with 


constant density cores (similar to globular clusters without core¬ 
collapse) and concluded that less-luminous galaxies with steeply 
rising slopes are those galaxies which host the most massive 
black holes when normalized to the mass of the prior constant 
density core. This leads to wrong predictions (e.g. by a factor 30 
for NGC 4889) because it underestimates the merger history of 
these galaxies. Instead, SMBH scouring leads in a natural way to 
the depletion of galactic centers ([Ebisuzaki et al.|1991]|Milosavl-| 


jevic & MetTitt|20OT| Merritt|20b6 1 and it predicts that the most 

massive black holes reside in galaxies with the most pronounced 
cores. [Ravindranath et al.| ( |2002| ) treated both scenarios (adia¬ 
batic and SMBH scouring) as competing models and argued that 
black hole merging is superior in explaining the observed param¬ 
eter relations. Today, adiabatic black hole growth models have 
disappeared mostly from the literatur^ 

Motivated by the wide distribution of the central cluster 
black hole growth rates (Section]^, we hope to rehabilitate adi¬ 
abatic SMBH growth models - but this time not as a competitor 
model to SMBH scouring. The relocation of mass i.e. adiabatic 
black hole growth (|Young|1980| Goodman & Binney|1984 


van 


der Marel|1999|l leads to a cusp regeneration process of the de 


pleted density centers after SMBH scouring. Of course, the ob¬ 
tained results depend on initial conditions and the accretion rate 
of hot cluster gas. In this sense, our findings are not free of sys- 
tematics, too, as we also have to make assumptions about ini¬ 
tial conditions which might turn out to be inappropriate. How¬ 
ever, initial conditions can be modified in future studies. The 
method itself allows to construct a black hole mass gain sensitive 
“calorimeter" which is based on the innermost surface brightness 
profile and which can be used to estimate (independently of e) 
the amount of swallowed mass in galaxy clusters. 


As already motivated in Section |23| we must assume that 
black hole growth (through hot gas) is powered by accretion 
flows which extend far outward, e.g. classical Bondi or giant ad- 
vection dominated flows. Accretion from a nearby molecular gas 
reservoir which is not replenished can not lead to adiabatic cusp 
regrowth owing to an unchanged potential. Star formation which 
can also increase the central surface brightness is neglected and 
therefore our results only apply to the density profile of old stars. 
Furthermore, we use a simplified merger history and assume that 
core profiles are created by the currently most accepted model, 
the SMBH binary evolution scenario. In the scenario of cosmic 
structure formation, several major merger events (4-5) are ex¬ 
pected to occur between z = 0 and z = 3 ( Conselice|2007 1 with 
a strong decline at z = 1. Our models start at this redshift and 
we assume that core-profiles were fully grown up to z = 1, so 
we can use current parameter relations for fixing the break ra¬ 
dius. We neglect the 1 major merger events since z = 1. Al¬ 
ternatively, we assume that alternating phases of core creation 
through merging and adiabatic cusp regrowth are dynamically 
equivalent. 


During our referee process Jingade et al. ( 20 15^ presented new adia¬ 
batic growth models by using classical Sersic bulges for the initial setup 
of galaxies. Their models cover the other extreme - the least luminous 
galaxies with nuclear star clusters. 


3.1. N-Body Setup 

3.1.1. Runs 

Nine representative SMBH growth models from Section are 
selected for our A-Body computations. On a logarithmic scale, 
they equally cover the whole parameter space including the ex¬ 
treme scenarios. The models are listed with their initial param¬ 
eters in Table The most extreme growth mode (Model 9) is 
computed by means of A-Body computations only from red- 
shift z = 1 up to z = 0.1. In this way, we enlarge the number 
of promising clusters hosting such massive SMBHs within our 
light-cone. 

Fully spin powered models with e = 1 (i.e. AM, = 0) are 
automatically computed in form of our reference models. They 
do not lead to adiabatic cusp regrowth since all binding energy is 
radiated away. Note that the conversion of given cosmic epochs 
into redshifts is done with the public available “Cosmological 
Calculator for the Flat Universe" from Nick Gnedin by using 
Go = 0.3036 and Hq - 68.14km s 'Mpc '. 


3.1.2. Initial galaxy profiles and parameters 


We use core- Sersic models ( [Graham et al.|2003) [Trujillo et al.[ 
2004|l for our initial galaxy setup since we assume that galactic 


cores are created within the first 5-6 Gyr, prior to the onset of 
AGN cluster feedback. The models are generated in equilibrium 
(including central black holes) using the procedure described in 
Brockamp et al. (2014[). In order to exclude any dynamical in¬ 


stabilities related to anisotropic models (i.e. tangetially biased 
within Rbieak and radially biased at large radii), we only use mod¬ 
els with an isotropic velocity distribution. Following [Graharn^ 
[ak] ( [2003j l and [Trujillo et al.[ ( |2004[ l, the six parametric Core- 
Sersic models have a light profile, I (r), (i.e. 2D density profile 
assuming a constant mass to light ratio) of the form: 


( 6 ) 



Here, b is a normalization parameter such that half of the total 
luminosity is produced within the projected half-light radius, Rg. 
The parameter, b, is obtained by numerically solving the integral 
equation AlO in Trujillo et al.[P004) l. The Sersic index, n, spec¬ 
ifies the light concentration of the outer profile beyond the break 
radius, Rbreak, (for short /?/,) where the profile transforms into a 
power-law with slope y. The parameter a controls the sharpness 
of the transition and Ib is the light intensity at Rh- Selected pa¬ 
rameters for our models are given in Table and are discussed 
below: 

By first fixing the initial SMBH mass to be M,[Mi)\ = 
1.5,5,10,15,20,40, Rb is then extracted from the parameter re¬ 


lation in Rush et al. (2013 1 . Although values as high as Rb - 
3.8 kpc seem far-fetched, there exist central cluster galaxies e.g. 
Abell 2261 and Abell 85 which have similar or even larger break- 
radii ( Postman et al.[[20T^ Lopez-Cruz et al. 2014| l. The effec¬ 
tive radii. Re, and dynamical masses, Mdyn, of our model galax¬ 
ies are based on real galaxies (NGC 1399, M 60, NGC 3842, 
NGC 4889, NGC 6166) with similar black holes ([Magorrian 
[enL][T9Ml [Lauer et al.[[2007bl [McConnell & Ma[[2013[ l. TTie 
SMBH masses {M,[Mq\ - 15,20) are within the one sigma 
confidence limits of NGC 4889. Hence, NGC 4889 is used as 
a representative for our highest mass models. For the uppermost 
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Table 1: Initial parameters of the model galaxies 


Model 

T[Gyr] 

init.M,[M 9 ] 

^dyn[^9] 

RJkpc] 

R&[kpc] 

AM, [Mg] 

Comment 

1 

8 

1.5 

400 

4 

0.15 

0.06 

OB-low, e = 0.42 

2 

8 

5 

770 

7 

0.4 

0.17 

OB-low, e - 0.2 

3 

8 

10 

1550 

18 

0.8 

0.38 

OB-low, e = 0.1 

4 

8 

10 

1550 

18 

0.8 

0.7 

OB-low, e = 0.057 

5 

8 

15 

1750 

25 

1.3 

4 

OB-medium, e = 0.1 

6 

8 

20 

1750 

25 

1.8 

6 

OB-high, e = 0.2 

7 

8 

20 

1750 

25 

1.8 

13 

OB-high, e = 0.1 

8 

8 

40 

2000 

68 

3.8 

40 

OB-medium, e - 0.01 

9 

6.7 

40 

2000 

68 

3.8 

120 

OB-high, e = 0.01 


Notes. List of parameters selected for our models/N-Body computations. “OB" stands for observation-based models (Section [2.2^ . Slope values 
y = 0.2, a = 4.25 and n = 4 (Equation]^ are used for all models. The huge initial SMBH masses in excess of 10'°Mo are motivated at the end 
of Sectio n|2.2| and in Section [23] While dynamical masses, Af^yn, and effective radii, R^, of our model galaxies are constrained by representative 
galaxies, the break-radius, R/,, rs taken from the scaling relation in |Rusli et al.| ( |201^ . SMBH masses during the M-Body computations are linearly 
increased up to AM. over the typical galaxy cluster lifetime ofT x S Gyr. The AM. values are selected to equally cover (in log-space) the whole 
parameter range of SMBH growth rates which are specified in Section]^ Model 9 is computed from z = 1 until z = 0.1. In this way the amount of 
potential candidates hosting exhaordinary massive SMBHs is increased over a sufficiently large light-cone. 


value {M,[Mg] = 40) we use NGC 6166. So far, it is the only 
source with a measured black hole mass based on stellar kine¬ 
matics which is comparable to that limit ( Magorrian et al.|1998| l. 
However, its mass has not yet been confirmed by other studies. 
For galaxy models 8 & 9 the bulge mass of NGC 6166 is rounded 
upwards to = 2000. 

We use the central slope parameter y - 0.2 which is assumed 
to be shaped by SMBH merging. This value is in accordance 
with several observed core-type galaxies in |Rusli et al.|p0T3] l 
and Dullo & Graham ( 2014| l and it results in a 3D density pro¬ 
file which declines faster than r This is an important con¬ 
straint for isotropic models as it guarantees the dynamical stabil¬ 


ity within the influence radius of the black hole (Tremaine et al. 
1994|l. The transition parameter, a - 4.25, is obtained from the 


arithmetic mean values from all galaxies in Rush et ak] P013| l 
and [Dullo & Graham (2014 1 except NGC 7768. This particular 
galaxy has a too sharp transition value, a —» oo, which would 
otherwise bias the arithmetic mean. All models and their param¬ 
eters are summarized in Table [T] 


3.1.3. A-Body computations 


The main computations are performed with the Muesli software 
for collision-less dynamics ( jBrockamp et al.| |2014| l on six 
nodes of the VLBI-computer cluster at the Max Planck Institute 
for Radio Astronomy. The calculations lasted more than a 
month. Muesli is a multi-purpose platform for the simulation 
of elliptical galaxies, their nuclear black holes and (not relevant 
for this paper) globular clusters systems. As a check we also 
performed a few computations by direct N-body simulations 
using a self-customized version of NBODY6 with N - 5 ■ 10^ 
particles ( |Aarseth||l999 2003 [ Nitadori & Aarseth||2012| l. This 
dual strategy allows to exclude potential error sources. These 
are either related to two-body relaxation processes (relevant 
for direct summation codes) where too small particle numbers 
would lead to a cusp regeneration process around a central 
black hole ( jBahcall & Wolf||1976| ) or potential fluctuations 
induced by discreteness noise in the self-consistent field (SCF) 
method ( jHernquist & Ostriker||199^ . The SCF formalism is 
implemented into the Muesli software. We carefully tested (by 
finding a converging solution) the required particle number, N, 
and the required order of the radial and angular base functions. 
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We use 10^ particles for the Muesli simulations. Furthermore, 
in order to decrease the computational cost and to reduce 
numerical errors (i) Model 1 (& 2) are not evolved for 850 (520) 
i.e. 8 Gyr but only for 50 A-Body time units (AM, remains 
unaffected). This strategy is applicable as we are interested in 
adiabatic processes only and want to suppress relaxation, (ii) 
During computations with Muesli, the overall potential/density 
is not updated in the Models 1-6 in order to reduce relaxation 
effects caused by discreteness noise. Model 7 is computed with 
(red-dashed line) and without (red-solid line) potential updates 
(PU). It is the model where back-reaction effects of the changing 
density profile (due to SMBH mass growth) start to become 
important. Therefore, Models 8 & 9 where AM, is of the order 
of the initial black hole mass itself are only computed with 
frequent updates of the potential. 


Despite dynamical processes like SMBH wandering, dynam¬ 
ical heating and relaxation driven cusp formation which oc¬ 
curs in NBODY6 computations (see e.g. Brockamp et al.|20TT i, 
the agreement with the Muesli simulations is nevertheless very 
good. The outcome of a NBODY6 computation is included in 
Figure]^ (Model 9, solid-blue line) for a comparison. To exclude 
(as much as possible) even unknown systematics we also com¬ 
pute all models without adiabatic black hole growth (i.e. with 
AM, = 0) and use them as a reference and for comparison. They 
trivially correspond to fully spin powered models, too. Finally, 
the amount of swallowed particles by the adiabatically growing 
black hole ranged from insignificant to small values when com¬ 
pared with the gas accretion rates and is therefore not discussed. 


3.2. Results 

Figure shows the dynamical response of the projected (2D) 
density profiles to adiabatically growing SMBHs obtained with 
the MUESLI software ( jBrockamp et al.|[2014| ). We use dimen¬ 
sionless model units, R,p 2 D^ for a better comparison. These mod¬ 
els can be scaled to physical units by Rphys = 1.34ReR and 

P2D {Rphys) = P2D (R) ^dyn/ {\.34Ref oc p (^phys), where p is the 
surface brightness profile using an appropriate mass to light ra¬ 
tio. All computations (red curves for adiabatic models and black 
curves for reference models without growing SMBH) contain 
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Fig. 3: Projected (2D) density profiles of all nine models. Red lines represent adiabatic SMBH growth models after T = 8 Gyr (Models 1-8) and T = 
6.7 Gyr (Model 9) whereas black lines show reference models without growing SMBHs (i.e. AM. = 0). The density slopes, y, are fitted within half 
the size of the initial break radii by using a simple power law approximation. Hence, they are influenced by the outer profile and we quote them with 
their sign. In ascending order of the computed reference Models 1-9 these are: - -0.26; -0.26; -0.3; -0.3; -0.28; -0.23; -0.25; -0.3; -0.3. 

The quantity Ay = y^f - yadiab ^ 0 results from a growing density cusp due to adiabatic SMBH growth. Evidently, accretion rates with AM. > 
M. (to) lead to a cusp regeneration process. Smaller accretion rates (Models 1-5) have no strong impact on the innermost 2D density slope. In 
panel 9 a pure Sersic n - A model is added for comparison (black dotted line). The solid-blue line corresponds to the NBODY6 computation of 
the same model. It is in excellent agreement with the solid-red line. The red-dashed line is similar but is scaled to a 1.5 times larger host galaxy 
mass, similar to the Phoenix BCG (McDonald et al. 2013^. 


10^ particles. The resolution scale in Figure ^is R^hys - 16 pc 
(Model 1) to /?phys = 270 pc (Model 9). iMe central density 
slope, y, is measured within half the initial break radius. Ob¬ 
tained uncertainties are based on assumed Poisson errors. Signif¬ 
icant cusp regrowth is defined over the quantity Ay = yref-Tadiab 


which must be larger than the sum of the individual slope errors. 
Evidently, noticeable cusp growth is only observed in Models 5- 
9 which means that the black hole accreted more than 25% of its 
initial mass. 
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Fig. 4: The (analytical) initial (i.e. AM.[Mg] = 0, black-dashed line) 
and final velocity distribution, - I - (cTg -i- of our most ex¬ 

treme Models 8 & 9. We use the 3D radius for the evaluation of the 
velocity dispersion components, Error values are obtained 

from the bootstrapping method. Particles on eccentric orbits are at¬ 
tracted more strongly towards the innermost center such that circular 
orbits, i.e. a tangetial biased region develops slightly outside but still 
very close to the center. This effect is relevant for the stellar dynamical 
mass measurement of the central black hole. 
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In Figure 1^ we plot the velocity distribution, ff - I - (cr^ -h 
(T^)/2(t^, of Model 8 & 9. In very good agreement with the out¬ 
comes of semi-analytical approaches (e.g. [Goodman & Binney| 
1984| Quinlan et al. 1995[ their figure 2) the isotropic velocity 


distribution around the pre-existing seed black hole, j3 (R) - 0, 
changes. Particles on eccentric orbits are dragged more effi¬ 
ciently towards the innermost center where the velocity distri¬ 
bution stays isotropic. Slightly outside, a tangential biased re¬ 
gion which is dominated by particles on more circular orbits 
forms. This has consequences for the stellar dynamical black 
hole mass determination method. Extreme cool-core clusters like 
the Phoenix- or RX J1347.5-1145 cluster with rapidly growing 
SMBH masses (unless e > 0.1) are rare within our light-cone. At 
huge cosmic distances, the innermost isotropic region (if it has 
not already been destroyed by SMBH merging) might be out of 
the resolution scale and what is left is the region where circular 
orbits dominate. The measurable radial velocities will be biased 
towards lower values and thus give rise to an underestimation of 
the true black hole mass if not properly taken into account. The 
effect is enhanced, if the initial model already starts with a tan¬ 
gentially biased core as expected from SMBH scouring models. 


4. Discussion 

General results, the method to unveil the most massive central 
cluster SMBHs and the outlook are part of this section. 

4. 1. How to unveil the most massive black holes and 
implications for scaling relations 

In §|^we compiled a set of SMBH growth rates which are based 
on the mechanical feedback energies of central cluster AGN. The 
adiabatic growth of hot cluster gas results in a cusp regrowth 
process if the ratio of mass gain to initial SMBH mass (at z = 1) 
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exceeds 25%. This ratio will be lowered if the initial slope is 
smaller than y = 0.2 or more generally if the total mass within 
the core is smaller than assumed for our computations. Neverthe¬ 
less, SMBH growth in ordinary elliptical galaxies/BCGs (repre¬ 
sented by our Models 1-4) should not lead to a significant adia¬ 
batic contraction. This may change in the most massive and X- 
ray luminous clusters with the highest feedback energies. Their 
central cluster black holes likely grow by several billion solar 
masses since redshift one unless most of the binding energy is 
radiated away (e 1). Assuming accretion efficiencies even be¬ 
low e ^ 0.1, sources with the highest feedback energies might 
even double their masses during the hot cluster gas accretion 
phase. Therefore, some nuclear black holes in the hottest and 
most X-ray luminous clusters might have accumulated (through 
hot gas accretion and S MB H coalescences) masses of the order 
of IO^'Mq (see Section|4^for some candidates). If the observa¬ 
tion based models (Equation]^ with small e are realized in some 
of these galaxy clusters, even the most prominent cores will be 
healed or even replaced by steep density cusps. In Model 8 the 
initial core-Sersic model is transformed back into a shape which 
resembles a classical Sersic profile whereas Model 9 develops a 
density cusp. This has consequences for the prediction power of 
the locally-calibrated M,-Rbreak, 4T,-Rcusp and M,-Ldef scaling 
relations. Black holes at the upper mass limit do not necessarily 
inhabit galaxies with the most pronounced light deficits/cores. It 
could also affects the My - Rcusp relation (|Lauer et al.||2007a i 
which links the absolute V-band magnitude of a galaxy to the 
size of the cusp radius obtained by a Nuker-profile fit. As light 
deficits are decreased by adiabatic contraction the most massive 
black holes might be located in BCGs with cored-density pro¬ 
files which are less pronounced than expected from the total lu¬ 
minosity of their host galaxy. Eirst observational support for our 
prediction might already have been found in [Savorgnan & Gra^ 
|ham| ( |2015] l. However, it is necessary to include additional, even 
more extreme cases to see whether their observed trend contin- 


Einally, one can also argue that BCGs with huge break radii 
and light deficits like MS0735 which are know to be in a rapid 
mass accretion process (unless e ^ 1) ( [McNamara et al.||2009| l 
must indeed be ultramassive, otherwise the ratio of AM, to the 
initial black hole mass would be too large and the core would 
have been replaced by a much steeper cusp. In this sense our 
computations also strengthen the case that galaxies with huge 
central cores witness the presence of ultramassive black holes. 
The ratio of black hole growth by merging (carving out shallow 
cores) to hot gas accretion (steepening core profiles) is crucial. 


4.2. Massive black hole candidates and a test 

Three promising galaxy clusters which might host very mas¬ 
sive black holes are RX J 1347.5-1145, SPT-CLJ2344-4243 (the 
Phoenix cluster) and Abell 2029. In this paragraph we will dis¬ 
cuss how they are related to our results and why they should 
be used for constraining hot gas accretion models. Additionally, 
they might help to uncover the limitations of various galaxy- 
SMBH scaling relation. 


- RX J1347.5-1145: This object at redshift z = 0.45 is one of 
the most X-ray luminous (Lx - b ■ lO'*^ erg s“' in the [2-10] 
keV band) and massive clusters known ( [Gitti & Schindler] 
[2004[ ). Depending on the accretion efficiency and assuming 
that its mechanical AGN feedback balances central gas cool¬ 
ing as supported by ( Churazov et al.|2002[ Gitti et ^.|2007] l, 
its nuclear black hole might have swallowed several IO^Mq 
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of hot gas (or even more if e < 0.1, § \2.2) . However, this 
is only one aspect. In their lensing analysis (|Kohlinger & 


Schmidt||2014l used 3D pseudo-Jaffe models to parameter¬ 


ize the two central galaxies. For the brightest cluster galaxy 
they deduced a velocity dispersion around cr « 600km s * (a 
value which is also supported by gas spectroscopy but which 
might not necessarily result from a virialized state |Sahu et al.| 
|1998) l and a 3D core radius in excess of a few kpc (it will be 
reduced in a 2D projection). If future observations confirm 
these values the case for an extremely massive black hole 
in the core of RX J1347.5-1145 would be strengthened. It 
would require a black hole merger rate at the upper end of 
the predictions of cosmological rich-cluster simulations 
porte & White|2015] ) to carve out a kpc sized core in a galaxy 
with a huge velocity dispersior0and to protect it against core 
healing by hot gas accretion. 

Phoenix cluster: This cluster is even more massive and X- 
ray luminous than RX J1347.5-1145 (McDonald et al. 2012| l. 
It could be one of the few objects in the visible universe 
where a model as extreme as our Model 9 might be real¬ 
ized and whose nuclear black hole might exceed IO^'Mq. 
Based on this scenario and our initial A^-Body setup we pre¬ 
dict that the old/red stellar density profile of the Phoenix 
BCG (red curves in Figure]^ should resemble a classic Ser- 
sic n = 4 model which becomes slightly steeper and more 
cusp-like within O.OR,,. Interestingly, the inner most data 
points (within ^ O.lRg) of the measured surface brightness 
profile in Figure 3 from [McDonald et al. ( |2013| l seem to lie 
above the Sersic n - A reference model. However, note that 
our A^-Body models are not fine-tuned to exactly match the 
properties of the Phoenix BCG nor its redshift. A more com¬ 
pact size of its BCG of Rg - 17 kpc ( [McDonald et al.|2013| l 
but otherwise identical initial parameters results in a larger 
initial mass/light deficit. This would necessitate even larger 
amounts of swallowed mass, AM., in order to steepen its 
central slope. Likewise, the applied bulge mass of its BGC 
has an influence too, as shown by the two red curves in Fig¬ 
ure [3] 

Abell 2029: The central host galaxy, ICllOl, is listed as 
one of the optically most luminous and largest cD-galaxy 
( [Dressler|[1979[ [Lin & Mohr[|2004| l. Therefore, a merger 
driven growth of its central black hole at the upper lim¬ 
its (M, ai [5 - 7] ■ 10 '‘’Mo) of cosmological predictions 
( [Laporte & Whit^[2015| l seems reasonable. Additionally, 
ICllOl is embedded within a massive and relaxed galaxy 
cluster ( [Lewis et al.|2003| l. If classical Bondi accretion with 
its quadratical SMBH mass dependency ( [Bondi[[T952[ our 
Section [23] l is realized, the nuclear SMBH should accrete 
at high levels. A precisely measured core size/light deficit 
might help to measure the ratio of SMBH scouring to adia¬ 
batic growth in one of the most extreme galaxies. 


Finally, Lopez-Cruz et al. (20141 report on the BCG in the 
Abell 85 cluster. So far it is the most extreme representative of a 
core-type galaxy with a me asured break radius (cusp r adius) of 
Rb - 20 kpc (Rg - 4.6 kpc) (Lopez-Cruz et al. 2014 ® Based on 
the extrapolation of the M, 


R\ 


break? 


- Rcusp and M. - L, 


^def 


^ The black hole radius of influence, Rinfl, decreases with growing ve¬ 
locity dispersion since Rinp oc Also note, that the simulated clusters 


injLaporte & WhiteH2015j are several times less massive (even) at z = 0 
than RX J13 47.5-1145 and SPT-C LJ2344-4243 at z ^ 0.5. 

; Bonfini et al.|i20l'5| take a different position by showing 


Note that 


that its stellar profile can be well fitted with a classical Sersic model 
without the need for a depleted core. 


relations , its central black hole m ight also have a mass around 


10"Mo (Lopez-Cruz et al. 


20141. 


It would be highly informative and important for constrain¬ 
ing the limitations of gas accretion- and SMBH merging pro¬ 
cesses, especially with respect to their cluster environments, if 
the nuclear black hole masses in these four objects could be mea¬ 
sured and compared with empirical scaling relations. In this way 
our prediction that the truly most massive black holes are not 
necessarily located in galaxies with the largest cores could be 
put on its first test. 


4.3. Outlook 

We now discuss future applications with the black hole mass 
gain sensitive ’calorimeter’, how its sensitivity can be increased 
and describe a non-feedback regulated SMBH growth channel 
which under right circumstances can be realized in the hottest 
and most massive galaxy clusters. 


4.3.1. Future applications 


A robust calorimeter requires detailed information about initial 
conditions. With a statistically meaningful knowledge about cen¬ 
tral surface brightness profiles of BCGs at the onset of hot gas 
accretion, in combination with deposited AGN feedback ener¬ 
gies, the calorimeter can be used to determine the significance 
of cold- (during the quasar phase) and hot gas accretion as well 
as SMBH merging. It also offers a tool to measure the accre¬ 
tion efficiency parameter. The calorimeter method is based on 
photometric- instead of spectral-imaging, therefore it is much 
easier to apply than the direct measurement of SMBH masses at 
these redshifts. A precisely calibrated SMBH mass gain sensitive 
calorimeter offers the potential to constrain the growth history of 
individual galaxies and their central black holes. 


We want to illustrate this idea on the basis of NGC 4889 and 
NGC 4874, the central dominating galaxies in the Coma clus¬ 
ter. While NGC 4889 is indeed the brightest member, it has a 
nearly flat density core with a slope y = 0.03 and a break ra¬ 
dius of Rb - 970 pc ( [Lauer et al.[[2007^ . Due to the flatness 
of its innermost surface profile, we do not expect, on the basis 
of the results presented in this paper, a significant phase of hot 
gas accretion thro ughout its life. Its centr al black hole of M, = 
(6 - 37) ■ IO^Mq ( McConnell et al.||2011 1 was likely grown dur¬ 
ing a massive quasar phase whereas subsequent SMBH merging 
activity has later shaped its core. NGC 4874 on the other hand 
has a much more diffuse appearance and is thus classified as a 
cD-type galaxy. If we apply our results to this galaxy, we find a 
different history. Its central black hole was not formed during a 
massive quasar phase but instead by the merger of several very 
massive progenitor galaxies. This gave rise to its extended halo 
as well as a very large break radius of Rb - 1730 pc ( [Lauer[ 
[et al.[[2007b| l formed by SMBH scouring. Its central cusp pro¬ 
file, y - 0.12, is steeper than that of NGC 4889. Assuming its 
initial profile was flatter in the past (e.g. comparable to that of 
NGC 4889), hot gas accretion was likely relevant in this galaxy. 
The recent infall of the galaxy group dominated by NGC 4889 
( [Gerhard et al.[[2007] l maybe quenched its activity through the 
disruption of the accretion flow. We can estimate the central 
SMBH mass in NGC 4874 from its break radius, the empiri¬ 
cal M, - Rb relation given in [Rush et ak] ( [2013[ ) and add the 
mass required to increase its central slope by Ay = 0 - 0.09 
(taken from our Model 7). However, the value AM.[M 9 ] 10 

relevant for the highest increase Ay = 0.09 must been taken as 
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an upper limit: The measured slope in NGC 4874 is even flatter 
than 7 = 0.2 (which is chosen for our galaxy computations) and 
the dynamical mass within its core is smaller (due to its more 
extended size and flatter slope) than assumed for our reference 
Model 7. Consequently, even moderate hot gas accretion with 
a value around AM,[Mg] ^ 1 might have shaped the central 
slope of NGC 4874. Its nuclear black hole might therefore be 
as massive as M, (20 - 30) ■ IO^Mq with a tendency towards 
the lower value, making it comparable to its neighbor. This mass 
estimate is also compatible with dynamical mass measurements 
performed in Magorrian et al. (|1998]l. 


much in excess of those typically observed (Dullo & Graham 
|2014| l. Nevertheless, we hope that future studies which aim to in¬ 
vestigate the influence of hot cluster gas accretion onto the most 
massive black holes in giant cool core galaxy clusters consider 
both scenarios: (i) Adiabatic cusp formation by hot gas accretion 
which is not restricted by the M, - cr relation (see e.g. |Hlavacek- | 
Larrondo et al. (2012|l and the results obtained in this paper) 


models. For the sake of simplicity we assumed that the grow¬ 
ing black hole has no influence on the accretion rate over typical 
cluster lifetimes of 8 Gyr. However, due to increasing evidence 
for the existence of black holes in excess of IO'^Mq and because 
of the generic dM,/df = yMl dependency of classical Bondi ac¬ 
cretion, this assumption might turn out to be wrong - at least for 
a few extreme objects. From a formal perspective, Bondi accre¬ 
tion should amplify itself and accelerate as: 


M. (f) = 


M. i t = 0) 

1 - yM, (f = 0) f ’ 


(7) 


4.3.2. Potential oscillations 

In our computations we assume that initial core formation pro¬ 
ceeds via well established SMBH merger scenarios and that the 
cores are fully grown at the onset of hot cluster gas accretion. 
Contrary to this assumption Martizzi et al. ( |2012| 12013) 1 show, 
that based on idealized computations (e.g. final SMBH masses 
following M, - cr), centrally flattened profiles can also be pro¬ 
duced by oscillations of the central potential induced by repet¬ 
itive AGN outbursts. However, their models predict core sizes 


However, care has to be taken. First of all (i), the gas reservoir is 
finite, (ii) the Bondi accretion parameters which are expressed as 
y are not necessarily constant over time subject to AGN feedback 
and (iii) by deriving the classical Bondi accretion rate (Equa¬ 
tion 0 one has to assume that the mass within the loss cone is 
small compared to the mass of the accretor. Nevertheless, by us¬ 
ing Equation |7] and the initial parameters from Section [Z3] with 
a black hole mass around 4 ■ 10'°Mo, the accretion rate is ex¬ 
pected to accelerate. After several billion years it would exceed 
a few percent of the Eddington rate where it likely changes to 
a radiatively efficient accretion state (|Gallo et al. 2003) Mac- 


and (ii) the destruction of density cusps by potential oscillations. 
Without the cold clumps which are used in |Martizzi et al.| ( |2013| l 
to mimic infalling galaxies and whose extended gas halos might 
be stripped more efficiently in hotter and more massive galaxies 
clusters, core sizes might be reduced due to less violent poten¬ 
tial fluctuations ( |Martizzi et al.|20l3] l. With smaller initial cores 
created by potential oscillations, adiabatic cusp formation by hot 
gas accretion might dominate at least in the most massive clus¬ 
ters hosting the most violently growing SMBHs. Core creation 
by potential fluctuations can also be implemented in future ver¬ 
sions of the Muesli software. In addition to the stellar/dark mat¬ 
ter (DM) profile a particle distribution which represents the den¬ 
sity profile of the ambient gas can be generated in equilibrium. 
These gas particles are not dynamically evolved forward in time 
but their masses are temporarily lowered within (e.g. bipolar) 
volumes to mimic cavities whose size depend on the equation 
of state of the gas and the jet power. Depending on additional 
parameters like the duty cycle of the AGN, longevity of cavities 
and the fraction of jet energy which goes into cavity production 
and shock generation, the stellar/DM particles (which are dy¬ 
namically evolved forward in time) will feel these potential os¬ 
cillation and react correspondingly. In addition to that the mass 
of the black hole is adiabatically increased as already studied in 
this paper. The implemented angular base functions (spherical 
harmonics) already have a topology which strongly resembles 
observed cavity systems. Eurthermore, the high resolution scale 
of Muesli allows to study the impact of potential oscillations 
deeply within the influence radius of the black hole, making it a 
reliable alternative to grid based codes. 


4.3.3. SMBH runaway growth 


In Section 2.3 we compared observed AGN feedback energies 
with those ones obtained from the prediction of classical Bondi 


carone|2003[ |Falcke et al.|[2004[ |Wu & Cao||2008| i. Depending 
on the broad-band spectrum of the quasai[]and the gas tempera¬ 
ture of the cluster it might either lead to inverse Compton cooling 
( Fabian & Crawford) 1990| l which decreases the gas temperature 
and would stimulate even higher accretion rates or to Comp¬ 
ton heating ( |Sazonov et al. 2004) 2005 i. Interestingly, power¬ 
ful quasars with luminousities around L ^ lO'*^ - 10^^ erg s“' 
(comparable to our model predictions) are reported to be lo¬ 
cated in the Phoenix ( |Ueda et al.|2013| l and CL 1821-1-643 cluster 
( [Walker et al.||2014 l. These are two of the strongest cool-core 
and most massive galaxy clusters. 


5. Summary 

This paper is about a new strategy to unveil the most massive 
black holes in the universe. SMBHs at the highest mass scale are 
important for constraining the limitations of empirical scaling 
relations which are useful for theoretical model building. With 
growing evidence that the most massive black holes are located 
at the centers of very massive galaxy clusters (e.g. Hlavacek- 
Larrondo et al.|20T2 1 we present a simple and effective method 
to pick out the most extreme candidates. 

Eor that purpose we first transformed (observed) AGN feed¬ 
back powers into central cluster SMBH growth rates (Section]^. 
Through a comparison with classical Bondi models we checked 
whether the initial black hole masses adopted for our A-Body 
computations are reasonable. We found a strong overlap between 
theoretically predicted and observationally observed feedback 
powers. This indicates that simple Bondi accretion (even when 
adapted for rotating gas atmospheres) can fuel the most powerful 
cluster AGN. 

In the main part of our paper (Section]^ we computed with 
the help of the Muesli software the dynamical response of the 
host galaxy to the growing SMBH. In this way we constructed 
a black hole mass gain sensitive calorimeter. In principle, a pre¬ 
cisely calibrated “calorimeter” offers a way to measure the con¬ 
tribution of cold gas accretion, SMBH merging and subsequent 
hot cluster gas accretion to the growth of the central cluster 
black hole. Eor our initial A-Body setup we used core-Sersic 


Thanks to our referee for pointing this out. 
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models to account for SMBH scouring as well as nine differ¬ 
ent SMBH growth scenarios, encompassing in logarithmically 
increasing steps the whole parameters space of compiled gas 
accretion models. We found that adiabatically driven cusp for¬ 
mation is significant in the most extreme growth models, but 
can be neglected if the ratio of accreted hot gas to initial black 
hole mass is around 25% or smaller. In this way our results also 
conhrm the expectation that galaxies with large central core¬ 
profiles (especially if they are located in galaxy clusters with 
huge deposited feedback energies and by assuming they are not 
mainly spin powered) witness the presence of extremely mas¬ 
sive black holes. Otherwise the core would have been replaced 
by a steeper central slope. However, in the most X-ray lumi¬ 
nous galaxy clusters with permanent AGN feedback powers as 
high as (Pagn ~ 10“^® erg s“') even the most pronounced ini¬ 
tial cores will be turned into cusps if accretion efficiencies are 
smaller than e = 0.1. We argued that the nuclear black holes in¬ 
side RX J 1347.5-1145, the Phoenix cluster and Abell 2029 might 
be among the most massive ones in the universe. Secure mass 
measurements of their central black holes would not only help 
to uncover limitations of scaling relation but could be used to 
constrain the limitations of gas accretion- and SMBH merging 
processes. They also might represent promising candidates for 
future applications with the Event Horizon Telescope ( |Doele- 
|man et al.|2009| l. The hunt for Gargantua is on. 
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